ON THE BEHAVIOR OF SOLUTIONS TO SCHRODINGER EQUATIONS 
WITH DIPOLE TYPE POTENTIALS NEAR THE SINGULARITY 
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Abstract. Asymptotics of solutions to Schrodinger equations with singular dipole-type poten- 
tials is investigated. We evaluate the exact behavior near the singularity of solutions to elliptic 
equations with potentials which are purely angular multiples of radial inverse-square functions. 
Both the linear and the semilinear (critical and subcritical) cases are considered. 
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In nonrelativistic molecular physics, the interaction between an electric charge and the dipolc 
moment D G M. N of a molecule is described by an inverse square potential with an anisotropic 
coupling strength. In particular the Schrodinger equation for the wave function of an electron 
interacting with a polar molecule (supposed to be point-like) can be written as 



where e and m denote respectively the charge and the mass of the electron and D is the dipolc 
moment of the molecule, see [TT] . 

We aim to describe the asymptotic behavior near the singularity of solutions to equations asso- 
ciated to dipole-type Schrodinger operators of the form 



in M. N , where N ^ 3, A = 2m ^j D l ; being |D| the magnitude of the dipole moment D, and d = D/|D 
denotes the orientation of D. A precise estimate of such a behavior is the first step in the analysis of 
fundamental properties of Schrodinger operators, such as positivity, essential self-adjointness, and 
spectral features. Such an analysis has been carried out in [7] for the case of Schrodinger operators 
with multipolar Hardy potentials based on the regularity results proved in [8] for degenerate elliptic 
equations. In a forthcoming paper, the authors will apply Theorem 11.11 to establish spectral 
properties of multi-singular dipole type operators. 
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We emphasize that, from the mathematical point of view, potentials of the form A fe^ have 
the same order of homogeneity as inverse square potentials and consequently share many features 
with them, such as invariance by scaling and Kelvin transform, as well as no inclusion in the 
Kato class. We mention that Schrodinger equations with Hardy-type singular potentials have been 
largely studied, see e.g. [TJ El El [10l \17\ [18] and references therein. 

More precisely, in this paper we deal with a more general class of Schrodinger operators includ- 
ing those with dipole-type potentials, namely with operators whose potentials are purely angular 
multiples of radial inverse-square potentials: 

_ _ a(x/\x\) 

in R N , where N > 3 and a 6 L°°(S Ar - 1 ). 

The problem of establishing the asymptotic behavior of solutions to elliptic equations near an 
isolated singular point has been studied by several authors in a variety of contexts, see e.g. [13] 
for Fuchsian type elliptic operators and [5] for Fuchsian type weighted operators. The asymptotics 
we derive in this work is not contained in the aforementioned papers, which prove the existence of 
the limit at the singularity of any quotient of two positive solutions in some linear and semilinear 
cases which, however, do not not include the perturbed linear case and the critical nonlinear case 
treated here. Moreover, besides proving the existence of such a limit, we also obtain a Cauchy 
type representation formula for it, see and (jHJ). 

We also quote 12J, where asymptotics at infinity is established for perturbed inverse square 
potentials and in some particular nonradial case. Holder continuity results for degenerate elliptic 
equations with singular weights and asymptotic analysis of the behavior of solutions near the pole 
are contained in [8]. 

As a natural setting to study the properties of operators C a , we introduce the functional space 
D 1,2 (R ) defined as the completion of C^°(R JV ) with respect to the Dirichlet norm 

v l/2 

Mlx>i.2( R iv) := ( / |Vu(x)| 2 dx 



In order to discuss the positivity properties of the Schrodinger operator C a in D 1,2 (M. N ), we consider 
the best constant in the associated Hardy-type inequality 



(1) Ajv(a) ■= sup 



\x\ 2 a(x/\x\) u 2 (x) dx 



u6 DMCR")\{0} / \S7 u ( x )\*dx 



By the classical Hardy inequality, An (a) ^ ny^P esssu Ps^-i a ; where esssupgw-i a denotes the 
essential supremum of a in S . We also notice that, in the dipole case, by rotation invariance, 
An(t3 • d) does not depend on d. 

It is easy to verify that the quadratic form associated to L a is positive definite in T) 1,2 (Mr) 
if and only if An (a) < 1. The relation between the value An {a) and the first eigenvalue of the 
angular component of the operator on the sphere § w_1 is discussed in section [2l see Lemmas 12.41 
and 12.51 More precisely, Lemma T2.5I ensures that the quadratic form associated to C a is positive 
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N -2 x2 



definite if and only if 

Hi > 

where fii = Hi( a iN) is the first eigenvalue of the operator — Agw-i — a(8) on S^ -1 (see Lemma 
I2.1[) . We denote by the associated positive L 2 -normalized eigenfunction and set 



N-2 //iV-2 x2 



(2) a = a(a,N):= — + j + ^ 

In the spirit of the well-known Riemann removable singularity theorem, we describe the behavior of 
solutions to linear Schrodinger equations with a dipole-type singularity localized in a neighborhood 
of 0. 

Theorem 1.1. Let C K be a bounded open set such that G and let q g L^ c (fl \ {0}) be 
such that q(x) = 0{\x\~^ 2 ~ e ^) as |x| — > for some e > 0. Assume that a 6 L°°(§ JV_1 ) satisfies 
Aat(o) < 1 and « 6 JT 1 (f2), u a.e. in Q, u ^ 0, weakly solves 

(3) - Attfg) - "(f/jfl) M (x) = g(x) u(x) mil, 

\ x \ 

i.e. 

f a(x/\x\) f 
Vu(x) • Vw(i) dx — / — : — 12 — u(x)«;(x) dx = / g , (x)u(x) , u;(x) dec, for all w £ H (Q). 

TTien t/ie function 



u{x) 



|x|^i(x/|x|) 



zs continuous in anc 



( 4 ) iim ■ ■ u{x ] . = / f R-°u{ne) + [ R 



2cr + N - 2 



q(s6)u(s9) ds 



l-R N-l+a 







/or all R>0 such that 5(0, i?) := {x G R w : |x| «C i?} C fi. 

We notice that ^ is actually a Cauchy's integral type formula for u. Moreover the term at 
the right hand side is independent of R. In the case in which the perturbation q is radial then an 
analogous formula holds also for changing sign solutions to ([3]) , see Remark 14.21 

If the perturbing potential q satisfies some proper summability condition, instead of the stronger 
control on the blow-up rate at the singularity required in Theorem ll.il a Brczis-Kato type argu- 
ment, see [2], allows us to derive an upper estimate on the behavior of solutions. For any q 1, 
we denote as L q (ip 2 ,VL) the weighted L 9 -space endowed with the norm 

IML^fi) == (J^ 2 * {x)\u{x)\ s dx\ , 

where 2* = is the critical Sobolev exponent and and <p denotes the weight function 
(5) tp(x) := |x| CT ViO£/M)- 
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The following Brezis-Kato type result holds. 

Theorem 1.2. Let SI C i w be a bounded open set such that € fi, a 6 L 00 ^^ -1 ) satisfying 
Ajv(a) < 1, and £ L s (ip 2 , SI) /or some s > AT/2. Then, for any SI' <s Q, i/iere exists a positive 
constant 

C = C(N, a, ||F|| dist(Sl', 90), diam Si) 

depending only on N, a, H^H/,*^*™, dist(Sl', <9S1), and diamSl, suc/i i/iai /or any weak H 1 (fl)- 
solution u of 



(6) 



Au{x) - u {x) = p r - 2 (x)V(x)u(x), in SI, 



i.e. for any u £ H (0) satisfying 
V«(s) • Vui(x) dx 



/n M 



^ ^ 2 - u(x)w(x) dx = I (p 2 2 (x)V{x)u(x)w(x) dx, 



n 



<C||n|| L2 . (n) . 

The Brezis-Kato procedure can be applied also to semilinear problems with at most critical 
growth, thus providing an upper bound for solutions and then reducing the semilinear problem to 
a linear one with enough control on the potential at the singularity to apply Theorem 11.11 and to 
recover the exact asymptotic behavior. 

Theorem 1.3. Let CI be a bounded open set containing 0, a £ L°°(S A '~ 1 ) satisfying A^(a) < 1, 
and f : SI x R — > M such that, for some positive constant C, 

f(x,u) 



^C(l + M 2-2 ) for a. e. {x,u) £ SI x 



Assume that u £ _ff 1 (Sl), u ^ a. e. in SI, u ^ 0, weakly solves 
(7) 



A«(a:)-^M u ( a; ) = /( a;)U ( a; )) j 



i.e. 



Vu(i) • Vu>(x) da; - 
Then the function 

is continuous in SI and 
u(x) 



— j — — w(x)u;(x) dx = / f(x,u(x))w(x)dx, for all w £ Hq (SI) . 
u(x) 



\x\ a ipi{xl\x\) 



(8) lim 



x|->0 |x| CT Vl 



S«-i 



r _<7 u(r0) + 



-2a-N+2 



2o- + 7V-2 

e 7V-l+cr 



o 2a + N-2 



f(s6,u(s6)) ds 

f(s9,u(s8)) ds Ui(0)dV(0), 



/or all r > such that 5(0, r) := {x G 



pJV 



;| r} C SI. 
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Notation. We list below some notation used throughout the paper. 

- B(a, r) denotes the ball {x € M. N : \x — a\ < r} in with center at a and radius r. 

- dV denotes the volume element on the sphere S^ -1 . 

- un denotes the volume of the unit sphere S^ -1 , i.e. u>n — fgN-i dV{8). 

- for any a G i 1 (S Ar_1 ), /gjv-i a {9) dV{6) denotes the mean of a on S^ -1 , i.e. 

/ a{0)dV(0) = — [ a(6)dV{8). 

- the symbol ess sup stands for essential supremum. 

2. Spectrum of the angular component 

Due to the structure of the dipole-type potential of equation (|3|) , a natural approach to describe 
the solutions seems to be the separation of variables. To employ such a technique, we need, as a 
starting point, the description of the spectrum of the angular part of dipole Schrodinger operators. 

Lemma 2.1. Let a £ L°°(S N ~ 1 ). Then the operator — A§n-i — a(0) on S^ -1 admits a diverging 
sequence of eigenvalues /Zi < //2 ^ ••• ^ /■**<•• ■ the first of which has the following properties: 

(i) [i\ is simple; 

(ii) /ii can be characterized as 

. |V SW -iV^)| 2 dV{9) - a(0)^(0)dV(0) 

(> Ml ^h^mo} J SN ^^(0)dV(0) 

(iii) /ii is attained by a C 1 positive eigenfunction ipi such that mingjv-i ip\ > 0; 

(iv) if, for some k £ R, a{0) = n for a.e. £ S^ -1 , then [i\ = —k. 

(v) if a is not constant, then — esssupgw-i a < /ii < — f&N-i a{0)dV{0). 

Proof. We prove assertion (v), being (i), (ii), (iii), and (iv) quite standard. 
Since the function tp = 1 satisfies 



\V S N-iiP(0)\ 2 dV(0) - a(0)ip 2 (0)dV(0) = - / a(0)dV(0), 

we deduce that n\ ^ — f SN _ x a(0) dV(0). In order to prove the strict inequality, we argue by 
contradiction and assume that [i\ = — f SN _ 1 a(0) dV(0). Then tp\ = 1 attains the minimum 
value in ^ but, since a is not constant, it does not satisfy equation — AgN-i^i — a(0)ipi = 0, a 
contradiction. We can thereby conclude that /ii < —f SN _ 1 a(0)dV(0). 

From [TH1 Lemma 1.1], and the optimality of the best constant in Hardy's inequality, it 
follows that 

/ S iv-i \V S N^(0)\ 2 dV(0)-esssnp SN ^aJ SN ^^ 2 (0)dV(0) 

Irn \Vu(x)\ 2 dx - ((^f 1 ) 2 + esssupgiv-i a) j RN ^j^- dx 
= inf 5^ = — ess sup a, 

thus proving the left part of the inequality stated in (v). □ 
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(11) 



The asymptotic behavior of eigenvalues [ik as k — > +00 is described by Weyl's law, which is recalled 
in the theorem below. We refer to O [15] for a proof. 

Theorem 2.2 (Weyl's law). For a E L 00 ^ 1 ^^ 1 ), let {/J,k}k^i be the eigenvalues of the operator 
-A S N-i-a(6) onS^- 1 . Then 

(10) ^ k = C(N,a)k 2/{N - 1] (l + o(l)) as k^ +00, 

for some positive constant C(N,a) depending only on N and a. 

The following lemma provides an estimate of the L°°-norm of eigenfunctions of the operator 
— Agw-i —a(0) in terms of the corresponding eigenvalues. For classical results about L°° estimates 
of eigenfunctions of Schrodinger operators we refer to [TBI §4] and references therein. 

Lemma 2.3. For a E L 00 ^ 1 ^^ 1 ) and k E N \ {0}, let ip^ be a L 2 -normalized eigenfunction of the 
Schrodinger operator — A§jv-i — a(0) on the sphere associated to the k-th eigenvalue fik, i.e. 

(-A SN -!M0) - o(9) M9) = Mfe M9), ™ S N ~\ 
\f SN - 1 \M9)\ 2 dV(9) = l. 
Then, there exists a constant C\ depending only on N and a such that 

|^)KCiKI L(A, ~ 1)/4J+ \ 

where |_-J denotes the floor function, i.e. lx\ := niin{j G Z : j ^ x}. 

PROOF. Using classical elliptic regularity theory and bootstrap methods, we can easily prove that 
for any j E N there exists a constant C(N,j), depending only on j and N, such that 

H^fcll , 2(jv-i) < C(N, j)(nk + \\a\\ L ^(s«-^y ■ 

W (Jv-i)-4(j-i) (S^- 1 ) 

Choosing j — [^-f^-\ + 1, by Sobolev's inclusions we deduce that 

M /2, (JV ^_-J J '_ 1) (gAT-1) ^ C a,a( S N-lj ^ L°°(S N - 1 ), 

for any < a < 2(1 — + L"^t^"J)i thus implying the required estimate. □ 

Arguing as in the proof of [181 Lemma 1.1], we can deduce the following characterization of 
Ajv(a). 

Lemma 2.4. For a E L 00 ^- 1 ), let A N (a) be defined in (Qp. Then 
(12) A N (a)- J sN -.a m \0)dV(e) 



^m(^)\ { o} J sn i iv^-^WI 2 dV(9) + J SN ^ ^(8) dV(6) 

We notice that the supremum in (fT2]) is achieved due to the compactness of the embedding 
i/ 1 (S JV_1 ) <^-> L 2 (S Ar ~ 1 ). As a direct consequence of the above lemma, it is possible to compare 
Ajv(a) with the best constant in Hardy's inequality 

4 _ J RN u 2 (x)\x\~ 2 dx 

(N-2) 2 ue x>i.!'^)\{0} Irn \Vu(x)\ 2 dx 

Indeed, if a is not constant, there holds 

4 f 4 
— — — T + a{6)dV{9) < A N (a) < ess sup a— — — , 

(N - 2y j S n-i gjv-i (N - 2y 
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whereas, if a(9) = k for a.c. 6 e E> N 1 and for some k£I, then A N (a) = 4k/ (N - 2) 2 . 
Let us consider the quadratic form associated to the Schrodinger operator C a , i.e. 

Q a (u):= [ \\7u(x)\ 2 dx- [ "fr/iyfr) dx _ 



The problem of positivity of Q a is solved in the following lemma. 
Lemma 2.5. Let a G L 00 ^ 1 ^^ 1 ). The following conditions are equivalent: 

i) Qa is positive definite, i.e. inf — , Sf 1 -. \.„ , > 0: 

uepi.2(K«)\ {0 }J KJV |Vu(a;)| 2 da; 

ii) A N {a) < 1; 

iii) [i\ > — ( N ^ 2 ) 2 where is defined in (0). 

Proof. The equivalence between i) and ii) follows from the definition of Ajv(a), see ([1]). On the 
other hand, [18\ Proposition 1.3 and Lemma 1.1] ensure that i) is equivalent to iii). □ 

Remark 2.6. We notice that in the case of dipole potentials, namely if a(9) = A t^t, then, rewriting 
ilty) in spherical coordinates and exploiting the symmetry with respect to the dipole axis, An (a) 
can be characterized as 

■d\_ cos6w 2 (9)d9 



A N ( A ^7— r ') = A sup 

\ \X\ J ,„cfTi/f 



w£Hl(o,n) J* \w'(9)\ 2 + (jV ' 2) 4 (jV ' 4) (sing)- 2 w 2 (g) 



d,0 



In dimension N — 3, a Taylor's expansion of A at (A ^4) near A = can be found in |llj . 

There is no explicit formula for the values An (A . When N — 4, it can be expressed in terms 
of Mathieu's special functions. The results of a numerical approximation o/Ajv(A^p) performed 
with a finite difference method are listed in table QJ which highlights how the dipole Hardy-type 
constant detaches the classical Hardy constant more and more as the dimension grows. 



N 


(A*(i)r=^ 


[Mt^P 




N 


(A W (l))- X = ^li 




3 


0.25 


1.6398 




7 


6.25 


19.0569 


4 


1 


3.7891 




8 


9 


26.7407 


5 


2.25 


7.5831 




9 


12.25 


35.7231 


6 


4 


12.6713 




10 


16 


46.0044 



Table 1. Some numerical approximations of An (t^t) obtained by finite difference with 10000 steps. 
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3. A Brezis-Kato type lemma 



In this section, we follow the procedure developed by Brezis and Kato in [2] to control from 
above the behavior of solutions to Schrodinger equations with dipole type potentials, in order to 
prove Theorem 11.21 

Let us consider the weight tp introduced in |(SJ) and define the weighted H^-space H^{Vt) as the 
closure of C°°(fi) with respect to 



By the Caffarelli-Kohn-Nirenberg inequality (see [5] and [I]) and the definition of ip, it follows that, 
for any w £ £>J,' 2 (fi), 



for some positive constant Cjv,a depending only on N and a. 

Lemma 3.1. Let Q be a bounded open set containing and v £ Hi(f2) (~l L q (p 2 q > 1, be a 
weak solution to 




and the space I?^,* 2 (£7) as the closure of C^°(£l) with respect to 




(13) 




(14) 



div(ip 2 (x)\7v(x)) — ip 2 (x)V(x)v(x) 



in n, 



where V £ L s (ip 2 ,fl) for some s > Then, there exists a positive constant C = C(a,N) 

depending only on a and N , such that for any O' <s fi, v £ L^ 1 (p 2 , 57') and 




q-(2-2* ) 




where C(q) :— min {i, — j-j} and 



X 
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PROOF. Holder's inequality and (fT3|) yield for any w € P*' 2 (0) 

(16) J <p 2 *(x)\V(x)\w 2 (x)dx J tp 2 * (x)w 2 (x) dx + J tp^ (x)\V (x)\tp 2 (x)w 2 (x) dx 

\V(x)\^e q \V(x)\^£ q 

2 

< tqj <p 2 * (x)w 2 (x)dx + (J <p 2 ' (x)w r (x) dx^j 2 ( J Lp 2 * (x)\V(x)\% dx 

\V(x)\>l q 

o T {x)w 2 (x)dx + C N Jj y 2 {x)\Vw\ 2 dx\( J ip 2 * (x)\V{x)\^ da 

\V{x)\>l q 

By Holder's inequality and by the choice of £ q it follows that 

N 

(17) J cp 2 *(x)\V(x)\^ dx^ (^J cp r (x)\V(x)\ s dx\ 2 " ( J (p r (x)dx 
\v(x)\^e q " \v(x)\^e q 



[ I V ?\x)\V[.r)\~ t l.r) [ J l^f^) r" (-r)'l-r 

)v(x)\^e q q 



\V(x)\ 



and hence from (jT6j) we obtain that for any w € V^ 2 (Q) 

(18) J tp 2 * (x)\V(x)\w 2 (x)dx ^ £ g J tp r {x)w 2 {x) dx +mmi^,-^-^\(J p 2 {x)\Vw\ 2 dx 
Let rj be a nonnegative cut-off function such that 

2 



supp(77) g(], ?y = 1 on ft' , and |V?j(a;)| < — 



dist(fi',an) 

Set u n := min(n, |u|) e X>i< 2 (ft) and test dHJ) with v(v n ) q - 2 rf E V^ 2 (n). This leads to 

(q-2) f <p 2 (x) V 2 (x)\Vv n (x)\ 2 (v n (x)) q - 2 dx + f ip 2 (x)r) 2 (x){v n (x)) q - 2 \Vv(x)\ 2 dx 
Jn Jn 

= [ tp 2 *(x)V(x)r] 2 (x)v 2 (x){v n (x)) q - 2 dx-2 [ p 2 {x)r l {x)v{x){v n (x)) q ~ 2 Vv(x) ■ Vrj(x) dx. 
Jn Jn 

We use the elementary inequality 2ab ^ l/2a 2 + 4b 2 and obtain 

(19) (q-2) f p 2 {x)rf{x)\Vv n {x)\ 2 {v n (x)) q - 2 dx + \ f p 2 {x)if(x){v n {x)) q - 2 \Vv{x)\ 2 dx 
Jn 2 J n 

< / ip 2 * {x)V(x)rf{x)v 2 {x){v n {x)) q - 2 dx + 4 / tp 2 (x)\ Vrj(x)\ 2 v 2 (x)(v n (x)) q - 2 dx. 
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Furthermore, an explicit calculation gives 
(20) 



V(K)*- 1 x;7 7 )| 2 < (g + 4)(g 2) {v n ) q - 2 r) 2 \Vv n \ 2 + 2(v n )"- 2 \Vv\ 2 r, 2 



+ 2(v n ) q - 2 v 2 \Vrf + ±—{v n ) q \Vr, 
Letting C{q) := min {\, ^}, from (19]) and (gDJ) we get 

(21) C(q) [ ip 2 {x)\V({v n )i^ 1 vr 1 ){x)\ 2 dx ^ 2(2 + (7(g)) / ( /3 2 (a;)(« ,l (a;)) <J - 2 z;(a;) 2 |V?7(x)| 2 da; 



<2 



C%)^— - / (^ 2 (x)(t;"(x))«|V77(x)| 2 da;+ / <p r (x)V(x)r ] 2 (x)v 2 (x)(v n (x)) q ~ 2 dx. 



n 



2 

Estimate (fT8j) applied to 77(1;™) i _1 w gives 

(22) / <f 2 ' {xWix^^^ix^^vix)] 2 dx ^ /" ^ 2 ( a; )|V(77(t;")*~ 1 z;)(a;)| 2 dx 

+ 4, / ^ 2 *(x)(7J™(x))«- 2 7J 2 (x)77 2 (^)dT 
JO 

Using ([22]) to estimate the term with 1/ in |[2"T]). ([T5]) yields 



n 



+ 4C j v, a (2 + C( g )) , ^^^n^ji^^iy^ji^ 



<?(<?) j n 

+ CV, a (<? - 2) / ^(^^"(^I^Vrj^)! 2 dx. 
Jo 

Letting n — > 00 in the above inequality, (fT"5)) follows. □ 

Proof of Theorem 11.21 Let w be a weak H 1 (O)-solution to ©. It is easy to verify that 
ip(x) := \x\ (T ^i(x/\x\) £ iJ 1 (fl) satisfies (in a weak i/ 1 (f2)-sense and in a classical sense in fi\ {0}) 

. . . a(x/\x\) , . 
-A<p(x) - tp{x) = 0. 

Then v :— — <E -H"^(f2) turns out to be a weak solution to (fTl)) . Let i? > be such that 

ft' <E Q' + S(Q,2i?) <E O. 

Using Lemma O in Oi := fi' + B(Q,R(2 - n)) <s ft' + B(0,2R), r x = 1, with q = qi = 2*, we 
infer that i> G (tp 2 and the following estimate holds 

11-11 <-> 2 2 , < ^(diam^^f 8 4 + %±j> + i^)" ll^l^. 2 . 

\C{qi) (Rri) z {RnY C{qi)J ^ ' 
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Using again Lemma O in Q 2 := 0' + B(0,R(2 - n - r 2 )) (e fix, r 2 = |, with q = q 2 = (2*) 2 /2, 
we infer that v <G L~^~ (ip 2 , Q 2 ) and 

t(2-2") 



(2 .,3 s^O 2 (diam0) 92 



8 4 4( g2 + 2) 2£ 92 \« 



C(diamO) 



cr(2-2*) 



91 52 / O 



4( 9l + 2) 2^ \« 



\C( qi )(R ri f (J?n) 2 C( 91 ; 



/ 8 4 4( g2 + 2) 2^ 2 \« 



Setting, for any n G N, n ^ 1, 
9 



™ = §(y) . n « :=^' + S^0,i?(^2-^r^, and r n = ^, 
and using iteratively Lemma 13.11 we obtain that, for any n £ N, n ) 1, 
(23) ||«[|i«»+i( V 2»no ^ II w IIl«"+i( v 2 *o„) 



sc ||w|| L 2« (y 2« n) 
We notice that 



C (diam 



<t(2-2*) 



r — n / 



4(^ + 2) 2£ 9fc V* 



fc=i 



(24) 
where 



n 



( 8 



VC( ft ) (i?r fe ) 2 {Rr k y 
4(g fc + 2) 2^ \t 



C(q k 



^\C(q k ) (Rr k f (Rr k ) 2 C(q k ) 



exp 



fe=i ■ 



= — log I 



16fc 4 



4fc 4 ( (?fe + 2) 21 



ik 



q k b \R 2 C(q k ) R 2 C{q k )y 

and, for some constant C = C[N,a, \\V\\ L s( v 2*^, dist(fi', Oil)) > 0, 



b k ~ 2 



log 



C 



1 (1 



2 V 2 



as fc — > +oo. 



Hence J^n=i ^« converges to some positive sum depending only on | V\\ L s^2*q^, dist(f2', 90), AT, 
and a, hence 



lim 

n — >+oc 



C(diamO) 



<r(2-2*) 



n 



4(g fe + 2) 2£ 9fc V fc 



^VC( ft )(flr t ) 2 (i?r fc ) 2 C( ft ), 

is finite and depends only on N, a, ||V|| i »( v ,2*n), and dist(Q', 9fi). Hence, from (|23|) . we deduce 
that there exists a positive constant C depending only on N, a, \\V\\ La ^ v 2' n ^, dist(0', dQ), and 
diamf2, such that 

IMU^+i^'fi') < C\\v\\ L 2' (ip 2- :Qj) for all n£N. 

Letting n — > +oo we deduce that v is essentially bounded in SI' with respect to the measure ip 2 dx 
and 

IML-O^n') < c IML 2 *( V 2 *,Q) = C Il u lli 2 *(f2)' 
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where IMIl 00 ^ 2 *^') denotes the essential supremum of v with respect to the measure ip 2 dx. Since 
ip 2 dx is absolutely continuous with respect to the Lebesgue measure and viceversa, there holds 
IMIi-ft^'n') = IMU°°(n')> hcncc v e L °°(W) and 

IMU°°(n') < C\\u\\ L 2* (n) , 
thus completing the proof. □ 

If the potential V in equation © belongs to L N ^ 2 (p r ,Q) (but to L s ((p 2 *,Q) for no s > N/2), 
although we can no more derive an L°°- bound for u/(p, we can obtain for u/ip as high summability 
as we like. 

Theorem 3.2. Let Q be a bounded domain containing 0, a £ L°°(Ei N ^ 1 ) satisfying Ajv(a) < 1, 
and V £ L N / 2 (<p 2 ,$V) . Then, for any Q' <<= and for any weak H 1 ^) -solution w to @|, there 
holds * £ Li (if 2 *, Q,') for all 1 < q < +oo. 

PROOF. The proof follows closely the proofs of Theorem 11.21 and Lemma \'S. II However, since we 
only require V £ L N / 2 (ip 2 , fi), we have that for any q there exists i q such that 

2*/ MT// \,w , ^ ■ { C n\ 2C V,a 

ip {x)\V{x)\ * dx <minj — — 
\v(x)\^e q 

but we can no more estimate l q in terms of q, as we did in (|17j) thanks to the summability 
assumption V £ L s (ip 2 , f2) for some s > AT/2. Hence we still arrive at an estimate of type (f2U|) but 
we have no control on the product in (|24p as n — > +oo. □ 



4. Behavior of solutions at singularities 

The procedure followed in this section to prove Theorem 1 1 . 1 1 relies in comparison methods and 
separation of variables. Indeed we will evaluate the asymptotics of solutions to problem ([3]) by 
trapping them between functions which solve analogous problems with radial perturbing potentials. 
To this aim, the first step consists in deriving the asymptotic behavior of solutions to Schrodinger 
equations with a potential which is given by a radial perturbation of the dipole-type singular term. 
In this case, it is possible to expand the solution in Fourier series, thus separating the radial and 
angular variables, and to estimate the behavior of the Fourier coefficients in order to establish 
which of them is dominant near the singularity. 

Proposition 4.1. Let a £ L 00 ^- 1 ) be such that A N (a) < 1, R > 0, and u £ H 1 (B(Q, R)), u^O 
a.e. in B(0,R), u^O, be a weak H 1 -solution to 

~a{x/\x\) 



(25) - Au(x) 



MM) 



u(x) inB(0,R), 



where h £ LJ^ c (0, R) n L p (0, R) for some p > N/2. Then, for any r £ (0, R), there exists a positive 
constant C (depending on h, R, r, a, e, and u) such that 

^\x\ a ^ u{x) C\x\ a for allxe 5(0, r) \ {0}, 
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where a is defined in (0). Moreover, there exists a positive constant C (depending on u, h, N , and 
a) such that, for any 9 € S^ -1 , 

(26) lim u( P 9)p- a = CMS), 
and, for any r G (0, R), 

(27) G = J (r-°u(rr,) + jf ^ s _^°_ - h(s)u(s rf) ds 

- r- 2 - w + 2 j( 2g + jv _ 2 ft(*)t*0» 77) ds) ^ (77) fo) . 

Furthermore, for any r G (0, i?), t/iere exists a positive constant C (depending on h, R, r, a, e, but 
not on u) such that 

(28) u(p6)^C\\u\\ m{B {o,R))p a i for all < p < r. 

Proof. Let r G (0, R). We can assume, without loss of generality, that R > 1 and 7- = 1. Indeed, 
setting ?/j(a;) := u(rx), we notice that w € i/ 1 (i?(0, R/r)) and weakly solves 

"o(s/M) 



-Atu(a;) 



7«(x) in 5(0, i?/r), 



where := r 2 h(rp) satisfies h £ Lf^ c (0, R/r) fl L p (0, R/r). Hence, it is enough to prove the 
statement for R > 1 and r = 1, being the general case easily obtainable from scaling. 

Let R > 1, r = 1 and u E ^(5(0, i?)), u > a.e. in S(0, R), u ^ 0, be a weak solution of ([25]) . 
By standard regularity theory, u G C°(B(0, 1) \ B(0, s)) for any s G (0, 1). For any fc G N \ {0}, 
let ipk be a L 2 -normalized eigenfunction of the operator — A§jv-i — a(6) on the sphere associated 
to the fc-th eigenvalue pk, i-e. satisfying (fTl) . We can choose the functions tpk in such a way that 
they form an orthonormal basis of L 2 (S Ar_1 ), hence u can be expanded as 



(29) u(x)=u(p6) = J2<Pk(p)M0), 

fc=i 

where p = \x\ G (0, 1], 6» = x/\x\ G S w -\ and 

(30) MP) = / u(p0)^*(0)dV(0). 

J§N-1 

The Parseval identity yields 



/ |K(p0)| 2 ^(0) = y> fc (p)| 2 , for allO 



<P<1, 

fe=i 



and hence 



L 2 (B(0 



„1 / OO s OO „1 

n«iii»fBfo.D) = / p"- 1 (£iv*(p)i 2 / / r - 1 iv*(p)i a «ip 



Equations ([25| and (fTT|) imply that, for every k, 



<fk(P) + — — -<Pk{p) - ^fk(p) = h(p)(p k (p), in (0, 1). 
P P 
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A direct calculation shows that, for some c*, c§ £ M, 



(32) <p k {p) = p°* ( <* + 
where 



1 8- a i +l 



P u k k 



h(s)ip k (s) ds +p CT M 4 + 



1 s-^+i 



h(s)ip k (s) ds 



(33) (T 



iV-2 



and ct,. 



7V-2 UN-2 



For the sake of notation, we set 
(34) 



+ r 1 s - a i +1 

A k {p) = p a x / — -h(s)^ k (s)ds 



and 

so that 
(35) 



B k (p)=p°k [c k 2 + 



1 s~^+ l 



p u k - u k 



h(s)if k {s) ds 



Mp) = c k 1 p a " + A k (p) + B k (p). 



Without loss of generality, we can assume that 



N 



N 

1 i 1 . fN-2 



v(^) 2 + ^ 



so, setting e = 2 — < £ < |(cr^ + ^p), f° r every A:. From Holder's inequality and Lemma [2~3l 
it follows that 

(36) |^*(p)|=p'J 

c lP °i\p k \\^\ +1 <- x 



1 s-< +1 



p a k -°k 



-h(s)[ / u{s6)ip k {p)dV{Q) ds 



gjV-1 



u) N dp 



+ i , « - 1 W - 1 



iV-1 . JV-1 



+1 



CT fc -°k 



\Lp(B(0,1)) 
U1 N 



I u IIl 2 *(B(0,1)) 



— + ~ \h{s)\\u{s6)\dV(9) ds 



i , 

(_ g -+ + l-.^_Ll-^L r l) . 2 p. , 

gV Ofc^ 1 p 2*p-2*-p ,j s 



W + 2 + | N-l | 



^n n cmi— \ 



+i 



|£p(B(0,1)) 
Wat 



|Lp(B(0,1)) 
Wat 



Il 2 *(b(o,i)) 



IML 2 *(B(0,1)) 



2*P /■ _+ iV-2 i n JV\ "I 1 
p2*p-2'-p \ a b ~ + Z --p-> 

2 * p ( N 2 I a + W ^ T ) 



1_ 1 

P 2^ 



2^ / + , JV^2 _ A 



where wn — JgN-i dV{9). In particular 

(37) A fc (/5)=o(p CT *) a S/9 ^0+. 
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Moreover 
(38) 



s- a x +1 h(s)<p k (s)\ds 
Jo 



JV — 1 JV-1 



gJV-1 



a nF i+: V|/i( s )|| u ( fl e)|dv'(^ ) ds 



N + 2 



iV 

iV + 2 



|£p(B(0,1)) 
Wat 



w IIl 2 *(b(o,i)) 

i|/ifc| L^J + i W*» W o.i)) 



/ - , 1 N-l N-l \ 2* p 
5V CT fe "t" A p 2* '2*p-2*-p ds 



OJN 



|L 2 *(B(0,1)) 



L2*p-2* 



JV-2 



)] > 



due to inequality <r fc < — N 2 2 . As a consequence, 



(39) 



1 s-^ +1 



lim 



h(s)ifk(s) ds 



h(s)cp k (s)ds, 



is finite. Since u G L 2 (B(0, 1)), from ([37]), fl39j, and the fact that p^- (6>) £ L 2 (5(0, 1)), we 
conclude that there must be 

(40) 

hence 

(41) 



1 s-°*+ 1 



°k ~ a k 



Bk( P ) = P a * 



-h(s)ifk(s) ds. 



Since u 6 C°(B(0, 1) \ B(0, s)) for any s G (0, 1), it makes sense to evaluate ifik at p = 1 and, from 
(l32l) and EDI), we have that 



(42) 



= p fc (l) + 



1 s- CT k +1 



^fc ~ CT fc 



h(s)(p k (s) ds. 



From J32l), EH), and El, we deduce that 



1 S -^+! 



(43) <p k {fi)=ff" + 



CT. -cr fc 



h(s)tp k (s) ds + 



-h(s)ipk(s) ds. 



1 



h(s)ipk(s) ds 



P 

lo H -°k 

From above, (|38p . and standard elliptic estimates (which allow to estimate u outside the singularity 
in terms of its -H^-norm) we obtain that, for some positive constant c depending only on N , i?, a, 
and h, 



(44) 



|ci < c\p k 



"J 



lff 1 (B(0,.R)) 



1 + 



2*P 



2*p-2*-p 
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Arguing as in (|38|) . we find that 
(45) 

W + 2 | JV-l | , , 



\Bk{ P )\ ^ 



|LP(B(0,1)) 



a k ~ a k UJn 

From ([55]) . (151?)) , and (|4"5"|) , we can estimate </3fc as 



m IIl 2 *(B(0,1)) 



[a=Ffe( e -^-^)] 



Af-2M i-^-Tf 



(46) 
where 



I / M / I fcl <r + , "felMfel^-"^ - ^ 1 ll /l ll-L!»(S(0,l)) || || e -E^l 



at 



W + 2 

L0 N N 



N + 2 
U! N Ol 



Claim 1: there holds 
(47) 



1 1 

p 2* 



JV-2y L -|-2^ 



|A fe (p)+B fc (p)K|^^ 



+ / ||^IUp(s(o,i)) 2 



: > \ / 4- — \ p- 1 

= 1 \UJN{(T k ~ CT fc ) £ — 



+ Oik 



where 



u lli 2 *(B(0,l)) 



,7a- 



/ 11^11^(3(0,1)) V" ,-, £ -~^a Tl 



W(o^ - cr fe ) 



n 



1 (N-2\? 

in—) 



Pk 



-1, 



i.e. jfc is the unique integer number such that 



Pk - £■ 



Notice that e < i(cr^ + -^) implies j fc ^2. 



To prove the claim, we observe that, from (|34|) and (|46|) it follows that 



(48) lAfcGo)! ^ \c\\p^ 



+ a k \pk 



+1 



II^I|£p(B(0,1)) 



jV p / p-l Qfg 



|LP(S(0,1)) 



u IIl 2 *(b(o,i)) i / + -x 

,,- ||/l|Uf(S(0,l)) 1 



- fife) (e-^-)V 



^1 



Il 2 *(b(o,i)) 



|LP(S(0,1)) 
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In a similar way, from (f4lj) and (|46[) we deduce that 
(49) I^WKIcJl^ 11 * 11 ^ ^ 1 



-^)[ (e + .+ _ CT -)_£ T ] £ F- 
ll^l+lll II / II^IU»>(B(0,1)) 



Summing up (jig)) and ([15]) , we obtain 



^ + — \ p~ 1 



fc w fc 



|A fc (p) + B fc ( P )i<i^i^ + i%Mi 4x 

£ • ' - ~ 

and hence, from (f35|) . 

(50) | yfc (p)|<|^(l+ ""'T^l ^r) 



„ / II^I|lp(s(o,i)) x " 

u IIl 2 *(b(o,i))( ~~~7~F 



2p : 



\uj N (a+ -a k )J ( a + + N_2 _ 2e ^ V 



«fc 



lMfe| L 4 J ll u lli 2 *(B(o,i)) 1 ; — — 



\uj N {a k -<J k )J (J(N=2 



2p 



2z- 



W(¥) ' + 1*-*)- 



Using (|50p . we can improve our estimates of Ak(p) and Bk{p) thus obtaining 



IHU*>(B(0,1)) 1 . / II^I|lp(B(0,1)) 



fe " fc ; / e " 



i i fei +in ii ( W h \\Lp(B(a,i)) 
+ a fc |Mfc|L * J + « i3 . (B(01)) + v f- 



p-i . 



and 



II^I|lp(5(0,1)) 1 / l|fe|Up(B(0,l)) 



op — 1 



+ Ofc 



+1 



I u IIl 2 *(s(o,i)) 

p 3e 



( ll fe llLP(B(0,l)) 
- CTfc ) 



JV-2 



(W+Mfc + 3e P J{E=iy +lih -2e 
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Summing up we find that 

\A k {p) + B k {p)\ < \c\\p< 

+ a k \p k 



l|fe|Up(a(o,i)) 2 / ||fe|Up(B(o,i)) 



I || / IMUp(B(0,1)) 

I u IIl 2 *(b(o,i)) — T~F =T 



4p 



,3e- 



N-2 



{^f+Pk-ie 



[^f + p k -2e 



An iteration of the above argument (j k ~ 1) times easily leads to estimate Claim 1 is thereby 
proved. 

Claim 2: the function s i— » s _<T fc +1 /i(s)(,Cfe(s) belongs to i 1 (0, 1) and 

(51) lim p-"<*ip k {p) = cf + / - + -h(s)ip k (s) ds. 

p~>o+ J Q a k -a k 

Indeed, from (|47[) . (f3"5)) . (f4"4")l . and the choice of j k , it follows that 

(52) \<p k (j>)\ < ^I|w||hx( B (o,i))^ £ "^, 

for some positive constant dfc depending on k (and on a, i?, h N). We distinguish now two cases. 
If j k e < J (^f^-f + p~k - e, then from ![3"4"|). and we derive that 

|A fc (p) + Sfc(p)| < d' k ll«l|HHB(0,l))P W * +1)E "^» 

for some other positive constant d' k depending on k (and on a, R, h, N), and hence, by (|44|) and 
the choice of j k , 



(53) 



| Vfc (p)| < 4'||u|Ux( B (0,l))P fe+1)E - 



Estimate (|53")) and the choice of j^, imply that the function s i— > s +1 ft-(s)<y9 J t(s) belongs to 
i^O, 1). Moreover, from (J4TJ) and d%3j) it follows that 



(54) 



IflfcOOl < 4"lkll^(B(04))P° fc+2) ^^ = o(/A + ) as p 0+, 



and the claim is proved. If j k e — y (^j^) 2 + Mfc — £ j then from (|34|) . (|4Tj) . and (|52|) we derive that 

|Afc(p)+S fc (p)| 7fclMliri(B(o,i))P CT * Ilogpl^, 
for some other positive constant j k depending on k (and on a, R, h, N), and hence, 

+ p— i 

(55) Im/O I < 7klMl.ff 1 (i?(0,i))/= >o ' h |logp| " ■ 

Estimate ([55)) implies that the function s t— > s _a fc +1 /i(s)y?fc(s) belongs to i 1 (0, 1). Moreover, from 
(HH) and it follows that 

(56) |5 fc ( (0 )|< 7 ; , || U || H1(s(0;1)) p^(^ P S ^-^+ e ^- 1 |log S |d S ) V = o(//*) asp^0+, 
and claim 2 is proved also in this case. 
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Let us fix k such that 

(57) i'^;^)) 2 ii^wo,!)) <i 4>4e and ± >aU v k> - k . 

u N {a k -<J k ) £ — UN£ *? y l(N=2f + IJLk 3 2 

From (|47|) and (f5T|) . it follows that, for all k fc and for some positive constant C2 (depending 
only on N, h, and a), 



\A k (p) +13^)1^14^ ^2 (-) 

i—1 

+ O 2 ||u||jji( B ( ,l)) — r — |Mfe| L 4 J 7 



—J + Mfc 
which yields, for all fc ^ fc, 

(58) |A fc (p) + S fc ( P )| < i| c ^|p< + 6 fcP ^-/ 2 , 
where 

5fe = C2||u|lffl(B(0 ' 1)) U.K + -^)J Kl c^i)i- 

From ([TO]). ^ fe ~ fc 2 /^- 1 ) and j fc ~ feVC^- 1 ) as _> +00 . Hence we have that 

(59) |&feKC3||w||H 1 (B(o,i))exp(-C 4 fc 1 /^), 

for some positive constants C3 and C4 depending only on AT, h, and a. In view of (f35|) and (158)) . 
we deduce that 

~|c£|p<^ < |^(p)| + \b k \p a t' 2 , for all fc > fc, 

and consequently 

- 00 «i 00 „i 00 „i 

iEi c ii 2 / P 2 ° t+N - 1 d P <2Y i j o w k (p)\ 2 P N - 1 d P +2Y,\b k \ 2 1 ft + »- x dp. 

Hence, from (IHT1) and (f551). we obtain that 



00 I k |2 00 I l 1 2 

(6°) E l^H-F < 8||tt||i, (fl(0il)) + 8 £ T7~T~Z+ < +00. 



From (HI 



k-l 

(6i) u{p9) P -"t = p -°i ipM^e) + J2p tTt ~ at p~° t <Pk(p)M0) + Y,P~ at *">(p)M0)- 

k=2 k=k 

From (I5T1) we deduce that 



(62) hm p-'i pi(p)Vi(0) 



1 S -<T++1 



c r 1 - cr 1 



V'iW 
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and 
(63) 



fc-i 

lim y2 P °t-° : tp-°itp k (p)ip k (0) = o. 



From (|3"5|) , ([55]) , ([50|) , (|55|l . and (J37J), we deduce that there exists some positive constant C5 
depending only on AT, h, and a, such that, for all p £ (0, 1/2), 

(64) ^p-^l^Cp)!!^^)] < C^^r ||cfl^ + |6 fe |p^/ 2 Jl^lLC^-D/^+i 



fc=/c 
„fcl2 



00 



J2(N + 2a+)p 2 ^-^\ f i k \ 2 ^ N - 1 ^^+ 2 



k=k 
L(AT-1)/4J+1 



k=k 



which implies 
(65) 



Una f]p-°t<p k (p)i> k (e) = 0. 

k=k 



Collecting dHTJ), (62]), Q, and (65]), we finally obtain that 



(66) 



lim u(p9)p CTl = 

P-+0+ 



1 s -a++l 



Ci - °"i 



— ft.(s)t/5i(s) ds 



We notice that, in view of 1(52]). ([3TJ]) . and (gDJ), 



(67) 



o"i — o"i 



— h(s)ipi(s) ds 



1 Z" 1 s~^ +1 

< Pi(l) _c 2+/ -+ -h(s)ip 1 (s)ds 

Jo o{ - a 1 



■ h(s)u(s 77) g?s 



^1(77) dV(?7). 



The limit in (|26p follows now from and (|6"T)) in the case r = 1 and by a change of variable 
inside the integral for r ^ 1. Moreover estimates (|4"4"]) . ([5U]) . ([55]) , ([54"]) . ([5^]) . the definition of jfc, 
and (|64p . imply that, for some C§ > depending only on N , R, h, and a, 



(68) u(p9)p 



r 1 s — 



v u k - u k 



k>i 



^ C 6 ||m|| h i (b 



(S(0,fi)), 
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for all < p < 1/2. On the other hand, standard elliptic estimates in B(0, 1) \ B(0, 1/2) yield, for 
some C7 > depending only on N, R, h, and a, 

(69) u{pd)p- a i ^C 7 \\u\\ H1(mR)) , for all - < p < 1. 

Estimate ([28]) follows from ([68]) and (j69|) . 



From (jnnj) and the positivity of u, it follows easily that c\ + J Q s + 1 _ h(s)(fi(s) ds 



< min ih\ < (9) ^ max^i, for all 6 § 

gJV-l gJV-1 



1 "1 
7V-1 



> 0. Since 



and, by standard regularity theory, u 6 C°(S(0, 1) \ £?(0, s)) for any s S (0,1), the proof of 
Proposition 14.11 will be complete if we show that 



(70) c{+ — r -/i(s)v5i(s)ds > 0. 



In order to obtain ([70]) . we need to prove the following 
Claim 3: if k S N \ {0} and 



1 S - fc + + ! 



k 



h(s)<pk(s) ds = 0. 



then <fk(p) = for all p E (0, 1). Indeed, if c\ + J Q h(s)<fk{s) ds = 0, then 



(71) Vk {p) = -p°k / _ -h(s)<p k (s) ds + p"* / — -h(s)ip k (s)ds. 

Jo <? k - <y k Jo ° k - <J k 

From claim 2 and (|51[) . we know that there exists a constant t k depending on k (and on a, R, h, 
u, N) such that 

Wk{p)\ < 4p CT * ■ 

Using the above estimate in (|71|). we can improve such an estimate as 

\<Pk{p)\ *S «fcP fc ^ + rr — ■ 

Using the above estimate in (|71j) . we can obtain the following further improvement 

\UN(<Jk ~<7 k )eJ (2-1) — 
Arguing by induction, we can easily prove that, for all j G N, 



\Vk{p)\ < tkp" k 



+ ( nh\\ L v {B{0 ,x)) v (p s y 



\uj N {a+ - a k )e J (ji) 2 
and letting j — > +oo, we deduce that ip k (p) = for all p <E (0, 1). Claim 3 is thereby proved. 
We are now in position to prove (|70[) . Arguing by contradiction, let us assume that 



(72) c{+ / — -h{s)<p 1 (s)ds = 

Jo 0i - a 1 
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and let fco > 1 be the smallest index for which 



+ 1 



CT fco 



-h(s)tpk (s) ds ^ 0. 



-h(s)(fk(s) ds — for all k, then <pk = 



Such a fco exists in view of claim 3; indeed if c\ + J Q ^ 
for all k and it would be identically zero, thus giving rise to a contradiction. Moreover, from (1721) . 
we have that fco > 1, and, by claim 3, (fk = in (0, 1) for all 1 ^ k ^ feo — !■ Repeating the same 
arguments we used above to prove (|66p . it is now possible to show that 



(73) 



lim u(p8)p CTfc « 



E 



°k 



— h(s)tp k (s) ds 



MO), 



where m ko is the geometric multiplicity of the eigenvalue /_tfc . We notice that the sum at the 
right hand side is a nontrivial function in L 2 (§ Ar_1 ) which, being fco > 1, is orthogonal to the first 
positive eigenfunction ipi. Hence the right hand side of ([73"]) changes sign in S N ~ X . Therefore the 
limit in (|73p implies that u changes sign in a neighborhood of 0, which is in contradiction with 
the positivity assumption on u. Condition (|70p follows and the proof of Proposition 14.11 is now 
complete. □ 

Remark 4.2. If we let the assumption of positivity of u drop, following the proof of Proposition ^ ■ 1\ 
we can still prove a Cauchy's integral type formula for u. More precisely, if u € H 1 (B(0, R)) is 
a weak solution to (|25[) in £>(0, R) which changes sign in any neighborhood of 0, with a radial 
potential h G L^ c (0, R) !~\ L p (0, R) for some p > N/2, then, following the notation introduced 
in (|33p and letting fco > 1 be the smallest index for which 



N-l 



T + 



k o u(rrf) 



T+ 

fco 



a k 



0Y 



■ h(s)u(s 77) ds 



KM 



— T k 



h(s)u(s V )ds)^ ko (r ] )dV(r 1 )^0, 



for any £ S and r 6 (0, R) there holds 



lim u(p9)p "^o = 



E 



{/c: {ik—f^ko } 



u(rr)) 



^fc 

i-o-r 



— h(s)u(s rf) G?S 



°"fc 



— h(s)u(si]) ds )ipk{r]) dV(rj) 



M°)- 



Without the assumption of radial symmetry of the potential, it is still possible to evaluate the 
exact behavior near the singularity of the first Fourier coefficient cpi (see (p?5|) and (I3U1) ). 

Lemma 4.3. Let a G L 00 ^" 1 ) be such that A N (a) <1, R> 0, and u e ^(3(0, R)), u ^ a.e. 
in 5(0, R), u ^ 0, &e a weafc H 1 -solution to 

~a{x/\x\) 



-Au(x) = 



q(x) 



u(x) inB(0,R), 
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where q G L£ c (B(0, R) \ {0}) f~l LP(B(0, R)) for some p > f . TTien, /or any < r < R, 



(74) limp— u{ P e)^ 1 {6)dV{e)= [r-°u(r6)+ - q(s 9)u(s 6) ds 

/•r „JV-l+cr 
-2<r-jV+2 / * 



o 2cr + TV - 2 



q{s6)u(s9) ds dV(0) 



Proof. The proof follows the lines of the first part of the proof of Proposition 14.11 By scaling, 
it is sufficient to prove |74|) for r = 1. Let 



u{x) = u(p6) = y^y fc (p)^fc(g) and q{x)u{x) =q(pd)u(p6) = ^ ( k (p)ip k (9) 

k=l 

where p = \x\ G (0, 1], = x/\x\ G S^" 1 , 



fc=i fc=i 

3JV-1 



MP)= u( P 9)M0)dV(0), Ck(p)= q( P 0)u( P e)MO)dV(0), 

and i/'fe is an L 2 -normalized eigenfunction of the operator — Agw-i — a(0) on the sphere associated 
to the fc-th eigenvalue p k , i-e. satisfying (TTT|) . The first Fourier coefficient ip\ solves 

V'l(p) + — <p'i(fi) % ^0°) = Ci(p) ^ (0, 1). 
A direct calculation shows that, for some c\, c\ G K, 

Vi(p)=p a i [c\+ I — —ti(s)ds)+ P °i j^ c i + J ___ Cl(s)ds 



where a+ = + ^/ (^) 2 + Mi and erf = - J (^r 2 ) 2 + Mi- From Theorem Q and 

standard regularity theory, we deduce that u(x) ^ const \x\ a i in B(0, 1), hence, Holder's inequality 
yields 

(75) / |s-^ +1 Ci(s)|ds< / s- a i +1 ( ( \q(s9)\u(s9)tp 1 (9)dV(9))ds 

Jo Jo V A"- 1 / 



^ const J s^^fj s^\q(s6)\ dV(6)^J ds 
^ const [[g||i*(B(o,i)) (^J s(1 ds^j <c 



In a similar way, we obtain that 

(76) s^ s _ ai+1 Cl(s) e £i(o,i). 

CT i 

Since u G L 2 * (-8(0, 1)), erf < of, from (J75|), (2U, and the fact that p"^ tpi(9) L T (5(0,1)), we 
conclude that there must be 

f 1 s- a i +1 
c 2 = - — +Ci(s)ds, 
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hence 

(77) p-^i Pl (p)=cl + 



i s -a++i 



-<T{ +1 



0"x 



— Ci( s ) ds for any p G (0, 1) 



Notice that, by standard regularity theory, <p\ is continuous at p = 1, thus, letting p — > 1 in (|77|) . 
we obtain 



(78) c\ = ^i(l) - 

Arguing as in (|75|) . we obtain that 



-(Tj +1 



(79) 



<?! 



:CiW 



■ Ci( s ) ds. 



ds ^ const p 



Since p > f , (77]), ([78]), and JTSJ) imply that 



lim o ffl <pi(p) = 99i (1) 

p-»0+ 



1 g-crf + l 



-Ci(s)ds + 



1 S -<- 



o"i - fi 



— Ci(s) ds, 



and l|74p for r = 1 follows. The result in the case r ^ 1 can be easily obtained just by scaling. □ 

In order to extend the result of Proposition 14.11 to the case in which the potential is a non radial 
perturbation of the dipole-type singular term, we will construct a subsolution and a supersolution 
which solve equations of type (|25p and the behavior of which is consequently known in view of 
Proposition 14. II 

Lemma 4.4. Let a G L 00 ^^ -1 ) be such that A N (a) < 1, CeM, and e > 0. Then, for all 



(80) 



< r < 



i^L(l-A„(a))J , z/C>0, 



4C 
-CO, 



and for all 7 G H 1 ^ 2 (dB(0, r)), 7 0, 7 ^ 0, t/ie Dirichlet boundary value problem 

~a(x/\x\) 



(81) 



-Au(x) = 



J(0,r) 



7, 



u(x), in B(0, r), 
on dB{0,r) 



admits a unique weak solution u G H (B(0, r)). Moreover u is continuous and strictly positive in 
B(0,r) \ {0}, and there exists a positive constant C depending on a, C, e, N , and r, such that 

( 82 ) IM|,F/i(B(0,r)) < C"ll7ll_H-i/2(a_B(0,r))- 



In addition, if 7 G W 2 - 1 / k ' k {dB{Q,r)) for some k > N/2, then u G C°(B(0,r) \ {0}). 

PROOF. For a fixed r satisfying (|80f and 7 G H 1 / 2 (dB(0, r)), 7 ^ 0, 7 ^ 0, let t; be the unique 
_ff 1 (i?(0, r))-weak solution to the problem 

'v G ^(5(0, r)), 
-Aw = 0, mB(0,r), 
u = 7, on dB(0, r). 
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By classical trace embedding theorems, it follows that 

(83) ||«||ffi(u(o,r)) < const(iV,r)||7|| H i/2 (as(:0jr)) , 

for some positive constant const(JV, r) depending only on JV and r. Let us define the quadratic 
form a : H%(B(Q, r)) x H£(B(0, r)) -> K as 



B(0,r) 



Q(w,it) : 
and $ £ .ff -1 (.B(0,r)) as 



f 



Vw(i) • Vw(x) — — — (a(x/|x|) + C\x\ e ) w(x)u{x) 



dx, 



B(0,r) 



a(x/\x\ 



C 



|2-e 



By Hardy's inequality, it is easy to verify that 



(84) 



Q(u,u) > 



1 - Ajv(a) 



4C*+? 



(JV 



v(x)u(x) dx. 



\Vu(x)\ 2 dx. 



B(0,r) 



Since ([50)1 implies that [l — Ajv(a; — ( N _ 2 yz 
coercive. Furthermore, the function x > a(x /|x|)|x'~ ~ 



— '— — 1 > 0, we conclude that the bilinear bounded form Q is 



■C\s 



belongs to L N + 2 (B(0, r)), hence 



$ is a bounded linear functional on Hq (B(0, r)). From the Lax-Milgram lemma we deduce that 
there exists a unique w G Hh (B(Q, r)) such that Q(w,u) = h-i li) „i for all u G Hn (B(0, r)). 
In particular w weakly solves 



(85) 




a(x/|x|) + C\x\ s w(x) = 



a(x/\x\) 



C 



12-6 



v(x), in B(0, r), 
on dB{0,r). 



Testing the above equation with w and using (|84| . Poincare's and Holder's inequalities and (|83|) . 
we obtain that 

IHI-ffMBCo.r)) < c(a,C,£ ) iV,r)||u|[Hi( B (o, r )) < c'{a,C,e,N,r)\\j\\ H i/2 {dB{0tr)) , 

for some positive constants c(a, C, e, JV, r) and c'(a, C, e, JV, r) depending on a, C, e, JV, and r. It 
is now easy to verify that u := it) + v G H 1 (B(0, r)) satisfies (|82|) and is the unique weak solution 
to (f8Tj) . Moreover, testing (f8Tj) with — u~ := — max{— u, 0} and using (|84|) . we obtain that 

1- Ajy(a) - 



= Q( u -,u-) > 



(JV - 2) 2 



B(0,r) 



|Vu~(x)rdx, 



which, in view of (|80|) , implies that w _ = a.e. in 5(0, r), i.e. u ^ a.e in 5(0, r). The Strong 
Maximum Principle allows us to conclude that u > in B(Q,r) \ {0}, while standard regularity 
theory for elliptic equations ensures interior continuity of u outside the origin. 

If, in addition, we assume that 7 G W 2 ~ 1/k > k (dB(0, r)) for some k > N/2, then v G W 2 > k (B(0, r), 
and hence v G C ' Q (5(0, r)), hence, from elliptic regularity theory applied to (|85|) outside 0, we 
obtain that u G C°(B(0,r) \ {0}). □ 



Proof of Theorem [T7Q Let R > such that 73(0, R) C 0. Since q(x) = 0(|x|-( 2 - £ )) as 
|x| — > for some e > and g G i{J c (0 \ {0}), there exists a positive constant C such that 
-C\x\-^-^ < g(x) s$ C|x|"( 2 - e ) for a.e. x G B(0,R). Let us fix f = f(R,N,q,a,e), such that 
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< r < min 



■^^-(1 — Ajv(a))] 1 ^ e |. We notice that the Maximum Principle implies that 

u > in B(0, R)\{0}, whereas standard elliptic regularity theory yields u £ W 2,k (B(0, R)\B(0, s)) 
for all s £ (0, R) and some k > N/2, and, consequently, u is continuous in B(0, R) \ {0}. Hence the 



function j r := u \ dB ^ r ^ belongs to W 2 1 / fe,fe (9i?(0, r)) for some k > N/2 and is continuous and 
strictly positive on dB(0, r) for all < r ^ r. From Lemma l4~4l we deduce that, for any < r ^ f , 
there exist w r £ H 1 {B(Q, r)) and ■u r £ H 1 {B{Q, r)) continuous and strictly positive in -B(0, r) \ {0}, 
weakly satisfying 

~a{x/\x\) 



-Au r (x) 



-i~\dB(0,r) 



and 




7r, 

a(x/|x|) 



CM 



+ C1x|- 2 + e 



u r (x), in £?(0, r), 
on dB(0,r) 

u r (x), in -B(0, r), 



7 r , on 3B(0,r). 

From Proposition 14. 1[ there exist two constants A-i > A\ > (depending on f, N, q, a, e, and u) 
such that 

(86) A 1 |x| <T < u f (x) and u f (x) < A 2 |x| CT , for all x £ B(0, r/2) \ {0}. 

Furthermore, for all < r ^ f , u — u r satisfies 

a(x/|x|) 



(87) 

(u — u r ) 
while u — u r satisfies 



A(u - u r )(x) 

9B(0,r) 



= 0, 



-A(u — u r ){x) 



a(x/|x|) 



+ C\x\- 2+e 



(u — u r ) 



dB(0,r) 



0. 



(u-w r )(x)>0, in£(0,r), 
on 9.0(0, r), 

(u-u r )(x)^0, in 0(0, r), 
on 90(0, r). 



Testing (|87|) . respectively ([88]). with — (u — w r ) , respectively (u — u r ) + , and using (|84|) . we obtain 
that, for any < r ^ f, 

(89) u r (x) < u(x) < u r (x), for all x £ 0(0, r) \ {0}. 
In particular, from and (|89p. we deduce that 

(90) Ai|x| CT < u(x) sC A 2 \x\< 7 , forallx£0(O,r/2)\{O}. 

Estimate (|90[) and the continuity of u outside the origin imply that there exists a positive constant 
C (depending on q, R, Q, a, e, and u) such that 

(91) ^\x\ a sC u{x) < C\x\ a for all x £ 0(0, R) \ {0}. 
Let us now fix S = 6(N, a, e) > such that 

(5 < min {e, /(^f+ii} 
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and set 



f = nnn{f, + \(™^) }, 

with C given in (|9ip . The function u defined as 
(92) u(x) = \x\°- 5 ip l (x/\x\) 

belongs to if 1 (.8(0, R)) and, for all < r ^ f , satisfies 



a(x/|x|) 



A,7(.n - TV 2 " " ^ ^( 2 /(^) 2 + Mi ~ <?) M~ 2 "(g) ^C|x|- 2+s u(x), in 8(0,r), 



Hence, for all < r ^ f, u — u r satisfies 

~a(x/\x\) 



on dB(0,r). 



-A(u — u r )(x) 



(<*-Ur)L B rtl.^ > °> 



.t • d + C|a;|- 2 + e 



(u-M r )(a;) > 0, in 5(0, r), 
on dB(0,r). 



\dB(0,r) 

Testing the above equation with — (u — u r )~ and using (|84[) , we obtain that, for any < r ^ f, 

(93) u r (x) «; w(x), for all x G B(0, r) \ {0}. 

From Proposition 14. \\ for any < r ^ f, the functions 

w r (x) u r (x) 
£ l-> i — ; — — ; — ;— ^ and x i ► — - — — — — — 

|x| ff Vi(^/NI) |a;| CT ^a ) 

have limits as a; — > 0, which, accordingly with (|26H27|) and taking into account the continuity of 
functions u r and u r up to the boundary \x\ — r, can be computed as 

u r (x) 



Lr 



= lim 



M->0 \x\ a <ipi(x/\x\) 

r~ a u(rrj) — C 



+ Cr- 2 °- N+2 



2a + N - 2 



s +e u r (srj)ds 



r JV-l+er 

s- 2+e u r (sv)ds)Mv)dV{r ] ), 



2(7 + N -2 



and 



L r : = lim 



u r {x) 



1*1-0 |a;| ff Vi(a;/|a:|) 

r~ a u{ri]) + C 



gJV-l 



,1— ff 



-2+£ 



2cr + AT - 2 



it r (s 77) (is 



r g N-U 



2(j + AT - 2 



From ((89j) and (f91^ . it follows that 



L„ = r 



gJV-l 



u{rrj) r tj)i{rj) dV(rj) + o(l) as r — > 0. 
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From (|93|) , ([92]) , and the choice of 5, we obtain that 



L r = r 



u(rr))ipi(rj) dV (r)) + o(l) as r — > 0. 



Hence, from Lemma T4. 31 we conclude that, for any i? such that 5(0, i?) C fi, 



(94) lim L r = lim L r = 

r— >0 r— >0 



— JJ-Sc— -fV+2 



Br a u{Rr)) + 



■ q(s rj)u(s rf) ds 



2a + N - 2 

s A-l+cr 

2g + N _ 2 Q(sv)u(sv) ds )Mv) dV{rj). 



In view of 



there holds that, for any < r ^ r, 

M r (x) 



^ lim inf ■ 



it(a;) 



~ r |I™o la;]""! ViWM) ^ \x\ a ^ip!(x/\x\) 

u(x) ,, w r (a;) 
^ lim sup - — ; ; — - — — ^ lim : — : — — - — — 

| x |_>o |x| a ^iVi(a;/|a;|) |x|->o \x\ a »i ipi(x/\x\) 

Letting r — > 0, we complete the proof. 



□ 



5. Behavior of solutions to the semilinear problem 

The L q and L°° bounds of solutions to dipole-type linear Schrodinger equations with properly 
summable potentials, derived in Theorems 11.21 and 13.21 allow us to obtain in the semilinear case 
analogous estimates. 

Theorem 5.1. Let ft be a bounded domain containing 0, o 6 L 00 (S JV_1 ) such that A^v(a) < 1, 
and f : x R — > K such that, for some positive constant C , 

f(x,u) 



^C(l + \u\ 2 - 2 ) for a.e. {x,u) £ Q x R. 

Then, for any Q,' <g SI and for any weak 1 (il) - solution u of ([7]). there holds — £ L°°(Sl') 

Proof. Let SI' <s SI and u £ if 1 (SI) be a weak 7? 1 (^-solution to 0. We set 

f(x,u(x)) 



V(x) 



ip 2 * 2 {x)u(x) 



and notice that u £ L 2 (SI) yields 



ip 2 *(x)\V(x)\ N/2 dx < +c 



Hence Theorem [3J] implies that ^ G L% 2 ,J1') for all 1 q < +oo. Since 



^ 2 *(a;)|T/(a;)| s ^ const 1+ / if 2 ' ' (x) 
n \ Jn 



u[x) 



<p(x) 



(2*-2)s 



we obtain that V G L s ((p 2 , 51) for all s ^j- 2 -- The conclusion follows now from Theorem 1 1.21 □ 
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Proof of Theorem 11.31 From Theorem 15.11 it follows that q(x) — is such that 

q G L™ C (Q \ {0}) and q(x) = 0(|a;|~( 2 ~ e )) as |x| — > for some e > 0. Hence the conclusion follows 
from Theorem ll.il □ 
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